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In E. coll, accurate cell division depends upon the oscillation of Min proteins from pole to pole. We 
provide a model for the polar localization of MinD based only on diffusion, a delay for nucleotide 
exchange, and different rates of attachment to the bare membrane and the occupied membrane. 
We derive analytically the probability density, and correspondingly the length scale, for MinD 
attachment zones. Our simple analytical model illustrates the processes giving rise to the observed 
localization of cellular MinD zones. 
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Understanding how proteins are directed to specific lo- 
cations within the cell is one of the key challenges of cel- 
lular biology. At the basic level, targeting of proteins 
to subcellular locations is governed by physical processes. 
A striking example is the system of Min proteins, which 
functions as an internal spatial oscillator in E. coli, and 
is necessary for accurate cell division [1]. The proper- 
ties and interactions of the three Min proteins, MinC, 
D, and E, have been revealed by recent experiments [2- 
11]. MinD is an ATPase — a protein which binds and hy- 
drolyzes the nucleotide ATP to ADP. In its ATP-bound 
form, MinD associates with the inner membrane of the 
cell [7, 8], where it recruits cytoplasmic MinC [6] and 
MinE [7] onto the membrane. Once on the membrane, 
MinE activates hydrolysis of ATP by MinD which re- 
sults in MinD dissociating from the membrane. MinE 
and MinD together produce a spatial oscillator with a 
period of ~40 seconds [7-10]. In each oscillation period, 
the majority of MinD molecules accumulate at one end 
of the cell forming a "polar zone" . The MinD polar zone 
then shrinks toward the end of the cell and a new MinD 
polar zone forms at the opposite pole. MinE is observed 
to form a ring at the medial edge of the MinD polar zone 
[12]. The MinE ring moves along with the polar zone 
edge as the MinD polar zone shrinks toward the end of 
the cell. The dynamics of MinC follows that of MinD 
[13]. Complexes of MinC and MinD on the membrane 
block the formation of a ring of FtsZ protein [4] , a neces- 
sary first step in determining the site of cell division [5]. 
The spatial oscillations of MinD and MinC from pole to 
pole thus ensure that an FtsZ ring does not form near the 
poles, so that cell division can only occur near midcell. 
Recently, several numerical models have been proposed 
to explain the oscillatory behavior and mechanisms of 
protein targeting [14-19]. 

One of the key emergent properties of the Min oscil- 
lations is the length scale for formation of a new MinD 
attachment zone. In filamentous cells {i.e. cells that grow 
but do not divide), several MinD zones are observed in 



a striped pattern with a characteristic length separat- 
ing MinD attachment zones [8, 20]. In this Letter, we 
study the processes giving rise to new MinD attachment 
zones in the cell. Using a simple model, we demonstrate 
analytically that the length separating MinD zones de- 
pends on (i) the cytoplasmic diffusion coefficient and the 
nucleotide-exchange rate of MinD, and (m) the rate of 
attachment of ATP-bound MinD to the membrane. 

Model. Our aim is to illustrate the physical processes 
of reaction and diffusion that give rise to preferred length 
scales in the Min-protein system. Our ID model for these 
processes is abstracted from the 3D numerical model of 
Ref. [18], which is based on measured properties of MinD 
and MinE and which reproduces the observed pattern 
in filamentous cells. Other models have also accounted 
for the oscillations of the Min proteins [14-16], including 
the striped pattern seen in filamentous cells. We begin 
by considering a fully formed "old" polar zone of ATP- 
bound MinD (MinD:ATP) in one half of the cell. MinE 
activates the ATPase activity of MinD, giving rise to a 
source of cytoplasmic MinD:ADP. Subsequently, there 
are two stages leading up to reattachment of MinD to 
the membrane: (z) diffusion of each MinD:ADP until nu- 
cleotide exchange transforms it into MinD:ATP, and (m) 
continued diffusion of the MinD:ATP until it attaches 
to the membrane (see inset to Fig. 1). The attachment 
rate of MinD:ATP increases with the local concentration 
of membrane-bound MinD; thus the old polar zone is 
much stickier than the bare membrane [18]. Our aim is 
to derive analytically the density of attachment p{x) of 
MinD:ATP outside the old polar zone, in order to obtain 
the length scale of the new MinD attachment zone. 

Filamentous cell, single source. We approximate a 
long, filamentous cell as an infinite ID line. The old 
polar zone of MinD occupies a: < 0, and bare membrane 
occupies X > (see inset to Fig. 1). To model the effect 
of the MinE ring at the edge of the MinD polar zone, 
we first consider that MinD:ADP dissociates from the 
membrane only at a: = 0. After stage (?) - diffusion of 
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MinD:ADP until nucleotide exchange - the probability 
density of each MinD:ATP is 



Pi{x) = / Pv 
Jo 



{x\t)Qi{t)dt, 



(1) 



where Qi(t) is the waiting-time distribution for single- 
step nucleotide exchange with average waiting time ti, 
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and P'p{x\t) is the distribution following diffusion in ID 
for time t with diffusion coefficient V, 
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Substituting Eqs. 2 and 3 into Eq. 1 yields the initial 
distribution of MinD:ATP in the cytoplasm 
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During stage {ii) — diffusion of MinD:ATP until mem- 
brane reattachment — we initially assume that the region 
of the old polar zone (x < 0) is infinitely sticky for 
MinD:ATP. Hence, if a MinD:ATP diffuses into a; < 0, it 
will immediately reattach and become part of the old po- 
lar zone. Thus only those MinD:ATP molecules formed 
at a; > which never cross the origin in the course of 
cytoplasmic diffusion will attach to form the new zone. 
Given an initial position a;o > of a MinD:ATP, the 
probability that it diffuses to some x > after time t is 
Pvix-xo\t). We define the probability that a MinD:ATP 
diffuses to some a; > after time t without crossing a:; = 
by P2{x,xo\t). The probability distribution P2{x,xo\t) 
can be calculated using the continuum version of the Re- 
flection Principle [21], 

P2{x,Xo\t) = Pv{x - Xo\t) ~ Pv{x + XQ\t). (5) 

The attachment density p(x) for the bare membrane 
at X > can now be obtained from 



p{x) 



dxo / dtP2{x,XQ\t)Pi{xo)Q2{t), (6) 
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where Q2{t) — (I/T2) exp(— t/r2) is the distribution of 
waiting times for membrane attachment with average 
waiting time T2. Integrating over t and xq in Eq. 6, we 
obtain (for ti ^ T2) 
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where the total probability of attachment to the bare 
membrane at a; > is 



Ptot 



(8) 
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FIG. 1: (Color online) Normalized attachment density 
p{x)/ptot (Eq. 7) for different average attachment waiting 
times T2 from 0.1s to 10s, with the diffusion coefficient 
V — 2.5/im^/s and average nucleotide-exchange waiting time 
Tl — Is taken from Ref. [18]. Inset: cartoon of attachment 
processes. The solid blue line shows a MinD:ADP which 
diffuses to the position xo before undergoing nucleotide ex- 
change. The resulting MinDiATP then continues to diffuse 
until it attaches to the bare membrane at x, without ever 
crossing the edge of the old polar zone at a:: = 0. In contrast, 
the dashed blue line shows a MinD:ADP which undergoes 
nucleotide exchange at the position xg, where it immediately 
reattaches to the membrane in the old polar zone. 



In the limit T2 — > 0, where the membrane is everywhere 
perfectly sticky for MinD:ATP, Eq. 8 gives ptot = 1/2, 
as required by symmetry. In the opposite limit, T2 ^ ti, 
where sticking of MinD:ATP to the bare membrane at 
a; > is slow, only a small fraction, ptot = -\/ri / 4t2 , 
attach to the bare membrane. 

We expect that polymerization of MinD:ATP on the 
membrane leads to autocatalytic membrane attachment 
that enhances the peak in the attachment density [18]. 
Thus, the relevant length scale to characterize the new 
MinD:ATP zone is the position of the maximmn in p{x), 
which occurs at 



\/VtiT2 / T2 
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In Fig. 1, we plot the normalized MinDiATP attach- 
ment density p{x)/ptot from Eq. 7. At large distances, 
the profile decays exponentially as expected for a diffu- 
sive process (see Eq. 4). However, for a simple diffusive 
process we expect the profile to have a maximum at the 
source position, i.e., at x = 0. In contrast, we find that 
the probability of attachment is zero at the edge of the old 
polar zone at x = 0. This "zero boundary condition" fol- 
lows from the infinite stickiness of the old polar zone for 
MinD:ATP. Any MinDiATP which forms near the old po- 
lar zone has a high probability of crossing into x < and 
immediately reattaching as part of the old polar zone. 
Only those MinDiATP which form sufficiently far from 
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the old polar zone are likely to reattach as part of the 
new MinD:ATP attachment zone. The two competing 
effects, the zero boundary condition at a; = and the ex- 
ponential decay due to diffusion as indicated in Eq. 4, set 
the length scale for the formation of the new MinD:ATP 
attachment zone. In the model from Ref. [18], the time 
scale for sticking to a bare membrane is T2 ~ 10s, giv- 
ing Zmax — 2.7/im. This length scale agrees qualitatively 
with half the center-to-center distance (~ 3/^m) between 
neighboring MinD attachment zones as observed in Ref. 
[8]. 

Filamentous cell, distributed source. In the preceding 
analysis, the polar-zone edge at a; = was taken to be the 
only source of MinD:ADP. However, in experiments some 
MinE is observed throughout the MinD polar zone [11, 
12], suggesting that cytoplasmic MinD:ADP is released 
throughout the old polar zone as well. What effect, if 
any, does this have on the length scale for formation of 
the new MinD attachment zone? Instead of assuming a 
single source for MinD:ADP at a; 0, we now consider a 
source with distribution w{xs), <0. This distributed 
source of MinD:ADP modifies Eq. 1: 
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where Xs describes the position of the MinD:ADP when 
it leaves the membrane. Integrating over i, we find for 
a; > 0, 



Pi(.x) 



dXsW(Xs 



(11) 

Comparing to Eq. 4, we see that a distributed source of 
MinD: ATP simply reduces Pi{x) by the constant factor 
in square brackets. Therefore, the attachment density is 
the same as in Eq. 7, with ptot reduced by the constant 
factor in square brackets. This leaves aimax unchanged, 
so the length scale for formation of the new MinD at- 
tachment zone is identical. 

Finite polar-zone attachment probability. How is the 
reattachment profile altered by relaxing the approxi- 
mation that the old polar zone is infinitely sticky for 
MinD:ATP? To treat the case where the mean attach- 
ment time to the old polar zone is ra > 0, we introduce 
the quantities pc{x,t) and pm{x,t), respectively denot- 
ing the densities of cytoplasmic- and membrane-bound 
MinD:ATP as functions of time. Our previous expres- 
sion for the membrane attachment density p(x) (Eq. 7) 
is then equivalent to pm{x,t oo) for the case — 
and X > 0. 

The time-dependent densities of MinD in the cyto- 
plasm and on the membrane are described by the fol- 
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FIG. 2: (Color online) Attachment density p{x) (Eq. 16) 
for different values of ra, the average attachment time of 
MinD:ATP in the old polar zone, with T) = 2.5/im^/s, av- 
erage nucleotide exchange waiting time t\ — Is, and average 
bare-membrane attachment time T2 — Is. Note the shift of 
the maximum toward a; = as ra r2. 



lowing equations 
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We are interested in the final density of MinD on the 
membrane Pm{x, t oo) for a; > 0, which, from inte- 
grating Eq. 13, is given by ^ Pc{x,t)dt = /(x)/t2. 
Similarly, we define Pm{x, t—too) = f{x)/T^ for a; < 0. 
Then the unknown function ,f{x) satisfies 
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Note that Pc{x, 0) is the initial distribution of MinD:ATP 
in the cytoplasm, which is given by Eq. 4. We can now 
readily solve Eq. 15 for f{x) in the regions a; < and 
a; > 0, subject to the boundary conditions that f{x) 
vanishes at ±oo and is continuous and diffcrentiable at 
a; = 0. 

The new MinD:ATP attachment profile for a: > with 
a finite binding rate 1 /t^ in the old polar zone is 



p{x) = 
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In Fig. 2, we plot p{x) for T2 = Is and four values of t^: 
Os, 0.1s, 0.4s, and Is. For T3 = Os, the density profile 
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FIG. 3: (Color) (a) The position of the maximum in the MinD-attachment profile, Smax (Eq. 17), as a function of T2 € [0s,3s] 
and ra £ [Os, Is], (b) The ratio of densities /(0)//(5;max). 



in Eq. 7 is reproduced. For T3 = T2 = Is, the profile is 
symmetric about x = with a maximum at x = 0. The 
intermediate values show a shift in the maximum toward 
X = as Ta approaches T2, and a finite density at a; = 0. 
The maximum occurs at 

""""^ y/T2 - yfrl \ 2T2y/TiT3 + 2x2X3 - TiTs + T1T2 / ' 

(17) 

which is less than Xmax (see Eq. 9) for all nonzero T3. An- 
other quantity of interest is the ratio of the MinD:ATP 
density at the polar-zone edge to the peak density, rj — 
f{^) I f{xraax)- This IS a mcasure of the spatial definition 
of the new polar zone. In Fig. 3, we plot (a) imax and (b) 
/(0)//(imax) for T2 € [Os, 3s] and T3 e [Os, Is], T3 < T2. 
Interestingly, except for very small r3 , imax is a function 
only of the ratio T2/T3. Note also that the ratio 77 is small 
only for very low values of T3. These results demonstrate 
that the asymmetry of binding probabilities to the old 
polar zone and to the bare membrane is crucial to estab- 
lish a new attachment zone well separated from the old 
one. 

In summary, we have studied the processes that give 
rise to the length scale observed in cellular MinD attach- 
ment zones. Our analysis indicates that this length scale 
is an emergent property of the system which depends 
primarily on (i) the rate of nucleotide exchange, {ii) the 
rates of membrane attachment, and [iii) the diffusion 
constant for MinD. It is important to note that these are 
biochemical properties of MinD which can be fine tuned 
by evolution. Our analysis suggests that these parame- 
ters are tuned to block minicelling at the cell poles and 
at the same time to permit cell division at the cell cen- 
ter. Finally, we note that the Min system serves as an 
example of how physical processes can produce subcellu- 
lar protein localization. While MinD is "targeted" to the 
cell poles, the mechanism is not based on recruitment by 
a biological target, but rather on a dynamical physical 
instability in the system. 
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